A New Signature of Dark Matter Annihilations: 
Gamma-Rays from Intermediate-Mass Black Holes 
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We study the prospects for detecting gamma-rays from Dark Matter (DM) annihilations in en- 
hancements of the DM density (mini-spikes) around intermediate-mass black holes with masses in 
the range 10 2 < M/Mq < 10 6 . Focusing on two different 1MBH formation scenarios, we show that, 
for typical values of mass and cross section of common DM candidates, mini-spikes, produced by 
the adiabatic growth of DM around pregalactic IMBHs, would be bright sources of gamma-rays, 
which could be easily detected with large field-of-view gamma-ray experiments such as GLAST, and 
further studied with smaller field-of-view, larger-area experiments like Air Cherenkov Telescopes 
CANGAROO, HESS, MAGIC and VERITAS. The detection of many gamma-ray sources not asso- 
ciated with a luminous component of the Local Group, and with identical cut-offs in their energy 
spectra at the mass of the DM particle, would provide a potential smoking-gun signature of DM 
annihilations and shed new light on the nature of intermediate and supermassive Black Holes. 
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I. INTRODUCTION 

Although many astrophysical and cosmological obser- 
vations provide convincing evidence for the existence of a 
"dark" component in the matter density of the Universe, 
the nature of this dark matter (DM) remains unkown. 
It is commonly assumed that DM is made of new, as 
yet undiscovered, particles, associated with theories be- 
yond the Standard Model of Particle Physics. Among the 
most widely studied candidates are the supersymmetric 
neutralino and candidates arising in theories with extra- 
dimensions, which appear difficult to constrain with di- 
rect searches (i.e. by looking for nuclear recoils due to 
DM particles scattering off nuclei) and whose prospects of 
discovery at future accelerators strongly depend on the 
details of the particle physics setup (for recent reviews 
see e.g. Refs. 0, 0)- Indirect searches via the detec- 
tion of annihilation radiation may provide an interesting 
alternative, but they are usually affected by large astro- 
physical and cosmological uncertainties. Furthermore, in 
many cases, the detection of an annihilation signal may 
be difficult to distinguish from less exotic astrophysical 
sources. An example of this is the case of the Galactic 
center, where such high-energy radiation has been re- 
cently observed by several different experiments, with- 
out providing any conclusive evidence for or against an 
interpretation in terms of DM annihilation products (see 
Refs. 0,0, IE 00 and references therein). 

Here we describe a scenario that may provide smoking- 
gun evidence for the annihilation of DM particles. If 
intermediate-mass black holes (IMBHs), with a mass 
ranging between 10 2 and 10 6 M Q (e.g. 0), exist in 
the Galaxy, their adiabatic growth would have modi- 
fied the DM distribution around them, leading to the 
formation of "mini-spikes" , that is, large, local enhance- 



ments of the DM density [Jj ■ The DM annihilation rate 
being proportional to the square of the number den- 
sity of DM particles, these mini-spikes would be bright 
gamma-ray sources, distributed in a roughly spherically- 
symmetric way about the galactic center, and well within 
the observational reach of the next-generation gamma- 
ray experiments. Their brightness and isotropy make 
them ideal targets of large field-of-view gamma-ray ex- 
periments such as GLAST [ljj. In case of a positive 
detection, Air Cherenkov Telescopes such as CANGA- 
ROO 11], HESS [13, MAGIC ,133 and VERITAS 
could extend the observations to higher energies and im- 
prove the angular resolution. We argue that the obser- 
vation of numerous (up to ~ 100) point-like gamma-ray 
sources with identical cut-offs in their energy spectra, at 
an energy equal to the mass of the DM particle, would 
provide smoking-gun evidence for DM particles. 

In this paper, we make predictions for the number of 
detectable black holes in two different IMBH formation 
scenarios. In the first scenario, IMBHs form in rare, over- 
dense regions at high redshift, z ~ 20, as remnants of 
Population III st ars, and have a characteristic mass-scale 
of a few 10 2 Mq_|15| (a similar scenario was investigated 
in Ref. |gL llti Tl7|). In this scenario, these black holes 
serve as the seeds for the growth supermassive black holes 
found in galactic spheriods 0|. In the second scenario, 
IMBHs form directly out of cold gas in early-forming ha- 
los, in a sense that will be specified below, and and are 
typified by a larger mass scale, of order 1O 5 M . We 
demonstrate that, with respect to Ref. 0, the latter sce- 
nario leads to qualitative differences in the mini-spike 
profiles with dramatic consequences for the detectability 
of gamma-ray fluxes. For both scenarios, we make de- 
tailed estimates of the population of IMBH in the Milky 
Way DM halo using a complete model of IMBH forma- 
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tion at high redshift, black hole mergers, and halo merger 
and evolution [l^. This allows us the unique ability to 
make a detailed study of the detectabilty of mini-spikes 
as gamma-ray sources. 

This paper is organised as follows. In Sec. II A, we 
review the evidence and formation scenarios for IMBHs. 
In Sec. II B, we describe the model that we employ to 
estimate the properties of the local IMBH population, 
and we present the main properties (radial profile, mass 
function etc.) of IMBH populations in Milky Way- like 
halos in Sec. II C. Sec. Ill is devoted to the calculation 
of the mini-spike profiles and Sec. IV to the DM annihi- 
lation fluxes. Sec. V contains our primary results on the 
observability of gamma-rays from the annihilation of DM 
around IMBHs. In Sec. VI, we discuss the implications 
of our results and draw our conclusions. We perform 
all of our calculations in the context of a standard, flat 
cosmological constant plut cold DM (ACDM) cosmology 
with J7m = 0.3, 17a = 0.7, h = 0.7 and a scale-invariant 
primordial power spectrum with a normalization set by 
cr 8 = 0.9. 



II. EVIDENCE AND PROPERTIES OF IMBHS 
A. The case for IMBHs 

In the last few years, observational and theoretical ev- 
idence has accumulated Q for the existence of compact 
objects, heavier than stellar black holes, but lighter than 
the so-called supermassive black holes (SMBHs) lying at 
the centers of galactic spheroids. We consider here the 
mass range 20 ^ Mimbh/M <, 10 6 , where the lower 
bound of the IMBHs mass corresponds to recent esti- 
mates of the maximum mass of the remnant of a massive 
stellar collapse |20j, and the upper limit roughly indi- 
cates the minimum mass of SMBHs, assumed to lie in 
the range 10 6 <S M SMBH / M <, 10 9 (see e.g. Ref. 
for a recent review). 

A hint of the existence of IMBHs is provided by the 
detection of bright, X-ray point sources, called ultra- 
luminous X-ray sources (ULXs) , that are a ppa rently not 
associated with active galactic nuclei pll l22l |23[ . Al- 
though many known X-ray sources are associated with 
neutron stars and black holes, this intepretation fails in 
the case of ULXs. ULXs would have to emit radiation 
far above the Eddington limit, if M 20M©, and their 
positions in their host galaxies are not compatible with 
masses M J> 10 6 M©, because dynamical friction would 
cause these objects to sink to the centers of their hosts 
on a timescale shorter than a Hubble time (e.g. Ref. 
Accretion by IMBHs has been advocated as a possible ex- 
planation jig. Another hint for the existence of IMBHs, 
although not conclusive, comes from stellar kinematics 
in globular clusters [24j|; the observed relation between 
the mass and the velocity dispersion in selected globular 
clusters may fall on the extrapolation of the analogous 
relation for SMBHs [HH 127||. 



From a theoretical point of view, a population of mas- 
sive seed black holes could help to explain the origin of 
SMBHs. In fact, observations of quasars at redshift z « 6 
in the Sloan Digital survey [28l I29I l3Cj suggest that 
SMBHs were already in place when the Universe was only 
~ 1 Gyr old, a circumstance that can be understood in 
terms of rapid growth starting from massive seeds (see 
e.g. Ref. |3j|). Furthermore, the growth of SMBHs 
through accretion and merging of heavy seeds may aid 
in the understanding some of the observed relationships 
between supermassive black hole masses and the proper- 
ties of their host galaxies and halos [3^.l33ll33.l35ll3^ . l37]| . 
Scenarios that seek to explain the properties of the ob- 
served supermassive black hole population generally re- 
sult in the prediction of a concomitant population of 
"wandering" IMBHs throughout massive DM halos and 
the intergalactic medium pj|,|3l,|3!|. However, despite 
their theoretical interest, it is difficult to obtain con- 
clusive evidence for the existence of IMBHs. A viable 
detection strategy could be the search for gravitational 
waves produced in the mergers of the IMBH popula- 
tion [llllillJIIllIllII, which may become possible 
with the advent of space-based interferometers such as 
LISA. 



B. IMBHs formation scenarios 

We focus here on two scenarios leading to the forma- 
tion of black holes at very different mass scales. In the 
first scenario (which we refer to as scenario A), black 
holes are remnants of the collapse of Population III (or 
"first" ) stars . Numerical simulations suggest that the 
first stars may form when primordial molecular clouds 
with ps 10 5 M cool by formation and distraction of H2 
into cold pockets at the centers of their DM halos, with 
typical densities of order 10 4 cm -3 and temperatures of 
order a few xlO 2 K 0,E3i an d become gravitationally 
unstable. 

Newtonian simulations suggest that the fate of 
Pop III stars is very different from the case of their 
metal-enriched, comparably less massive counterparts 
mentioned above. Zero metallicity Pop III stars 
with masses in the range M ~ 60 — 140 M and 
M <; 260 M Q collapse directly to black holes, stars 
with 140 < M/Mq < 260 M are completely disrupted 
due to the pulsation-pair-production instability, leav- 
ing behind no remnant, and again stars with masses 
M <; 260 Mt> collapse directly to black holes ( see a l so 
Refs. |48l l49i l50l l51|). The evolution timescale of these 
very massive stars is typically of order t* ~ 1 — 10 Myr. 
After this timescale, supernovae begin to explode, re- 
leasing energy and metals into the surrounding medium. 
In the standard picture of hierarchical structure forma- 
tion, the metal-enriched material will be collected at later 
times at the centers of more massive haloes, where new 
generations of stars will form. 

Interestingly, if a ~ 10 2 M black hole forms halos that 
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represent ~ 3c peaks of the smoothed density field, the 
resulting baryonic mass fraction in these objects would 
be comparable with the mass fraction in SMBHs [lflj . 
Additionally, such a scenario leads to the natural predic- 
tion of a population of "wandering" black holes in the 
halos of Milky Way-sized galaxies, with masses similar 
to their initial mass scale M ~ 10 2 M , as many of the 
relatively small halos [M ~ 10 7 M Q for 3er fluctuations at 
z = 18) that host early- forming black holes do not merge 
with the central gala xy, but orbit about the periphery of 
the halo [H| l38l l39j . We stress that black holes in this 
scenario may not necessarily form at the very centers of 
their initial host dark matter halos at high redshift, a 
circumstance that, as we shall see, may have important 
consequences on the detectability of IMBHs. 

To represent the predictions of this class of black hole 
formation scenario where black holes form at ~ 100 M Q 
from the remnants of the first stars, we use a model sim- 
ilar to that proposed by Madau and Rees and stud- 
ied in further detail by Islam et al. 38] and Volonteri et 
al. ■ Specifically, at z = 18, we populate halos that 
constitute 3<r peaks in the smoothed primordial density 
field with seed black holes of initial mass 100 M . We 
evolve these halos using an analytic model of halo growth 
that is focused on making many statistical realizations of 
the growth of a Milky Way-sized halo. After populat- 
ing progenitor halos at high redshift with black holes as 
described above, these processes of halo growth and evo- 
lution are treated as described in detail in [5^, Ej| and 
0,13- We refer the reader to these references for details 
and tests of the halo evolution models. For the purposes 
of this study, we take the mass of the Milky Way halo 
to be Mmw — 10 121 h^ 1 M and perform 200 statistical 
realizations in order to ascertain the expected range of 
observable IMBHs. 

The second scenario that we consider (scenario B) is 
based on the proposal of Ref. [37j and it is representative 
of a class of models in which black holes originate from 
massive objects formed directly du ring the collapse of pri- 
mordial gas in early- forming halos |55ll5^.l5^ . l58ll5^ . l60| . 
In this class of models, the initial black holes are massive 
(~ 10 5 M Q ) and the growth of SMBH proceeds in such 
a way that both mergers and accretion play an impor- 
tant role. We use the model of Ref. [53 to represent the 
predictions of models that start SMBH growth from very 
massive seeds. The proposal of Ref. [33 is as follows. 
During the virialization and collapse of the first halos, 
gas cools, collapses, and forms pressure-supported disks 
at the centers of halos that are sufficiently massive to 
contain a relatively large amount of molecular hydrogen 
(molecular hydrogen is the primary gas coolant in halos 
in the relevant mass range, see [6l| for a review). In halos 
that are both massive enough that molecular hydrogen 
cooling is efficient and which do not experience any ma- 
jor mergers over a dynamical time, a protogalactic disk 
forms and can evolve uninterrupted. An effective viscos- 
ity due to local gravitational instabilities in the disk leads 
to an effective viscosity that transfers mass inward and 



angular momentum outward 62] until supernovae in the 
first generation of stars heat the disk and terminate this 
process [33 . By the time the process terminates (of order 
the lifetimes of Pop III stars, ~ 1 — 10 Myr), a baryonic 
mass of order ~ 10 M loses its angular momentum and 
is transferred to the center of the halo. Such an object 
may be briefly pressure-suppor ted, but it eventually col- 
lapses to form a black hole [471 |63j . 

The requirements that the early-forming host halo be 
massive enough to form an unstable disk and that the 
halo not experience a major merger imprints a typical 
mass scale for halos within which this process occurs of 
order ~ 1O 7 M . In this case the characteristic mass of 
the black hole forming in a halo of virial mass M v is given 
by 

M bh = 3 .8xlO 4 M (^)( ! 4) 3/2 

(l0 7 M o ) ("IB lOMyr (1) 

where we have assumed that a fraction / is the fraction 
of the total baryonic mass in the halo that has fallen 
into the disk, z is the redshift of formation, k is that 
fraction of the baryonic mass which loses its angular mo- 
mentum that remains in the remnant black hole, and t 
is the timescale for the evolution of the first generation 
of stars [33 • The distribution of black hole masses is a 
log-normal distribution with a mean given by the charac- 
teristic mass above and a standard deviation OM bh = 0.9. 
The spread is determined by the spread in total angular 
momentum exhibited by halos of fixed mass in cosmo- 
logical A-body simulations of DM halo formation 
Using the prescriptions of the Koushiappas et al. [371 ] 
model, we can again populate halos with black holes at 
high redshift and evolve them forward to determine the 
properties of satellite black holes in a statistically large 
sample of Milky Way-like halos at z = 0. This is pre- 
cisely what was done in Ref. [l9| in order to study the 
gravity wave background and we refer the reader to this 
work for further details. As with scenario A, we take the 
Milky Way halo to have mass M M w = 10 121 /i _1 M at 
z = and construct 200 realizations of wandering black 
hole populations in halos of this mass. 

C. Intermediate-Mass Black Holes in Milky 
Way-sized Halos 

In the previous two sections, we outlined models for the 
production of IMBHs in the early universe and evolution 
of IMBHs in their host halos in the context of the hierar- 
chical CDM model of structure formation. Of course, as 
halos merge to form larger systems that eventually grow 
to the size of the Milky Way, black holes merge, produc- 
ing supermassive, central black holes and perhaps a de- 
tectable gravity wave signal. These products have been 
the focus of most previous work regarding these mod- 
els 0, [53, HU [13 ■ Consequently, these studies focused 
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FIG. 1: Mass function of unmerged IMBHs in the scenario 
B, for a Milky Way Halo at z~0. The distribution is based on 
an average of 200 Monte Carlo realizations of a halo of virial 
mass M v = 10 121 /i _1 Mq, roughly the size of the halo of the 
Milky Way. 



much attention on the merging of black holes as halos and 
galaxies merge. On the contrary, we are most interested 
in those pristine black holes that are orbiting within the 
Milky Way halo and have not merged with other black 
holes because these unmerged black holes may still reflect 
the properties of the dark matter density enhancement 
in which they formed. 

In scenario A, the mass spectrum of unmerged black 
holes is a delta function as described in Section II A. 
The average number of unmerged black holes per Milky 
Way halo is -/Vbh,A — 1027 ± 84, where the errorbar de- 
notes the la scatter from halo-to-halo. In scenario B, 
the total number of unmerged black holes per Milky Way 
halo is 7Vbh,B - 101 ± 22. We show in Fig. the final 
mass spectrum (i.e. at redshift z = 0) of black holes in 
scenario B. As expected the distribution follow closely 
the initial mass spectrum, with a characteristic mass of 
order « 10 5 M Q . The only deviation is that the over- 
all distribution is slightly broadened by the fact that not 
all black holes form at the same redshift in halos of the 
same mass (see Eq. pQ and Refs. The radial 

distribution of unmerged black holes is less trivial, and it 
would be more difficult to derive directly from the mod- 
els of IMBH formation at high redshift. The distribution 
is essentially set by the energy and angular momentum 
distributions of merging objects in a ACDM cosmology 
and dynamical friction (e.g. |53|). Unlike dark matter 
substructures, which are generally absent from the inner 
parts of the host halos, because they tend to lose mass 
via tidal mass loss and heating, black holes and the sur- 



FIG. 2: Cumulative radial distribution of unmerged IMBHs 
in the scenario A (red) and B (black), for a Milky Way Halo 
at z=0. The mean and error are based on 200 Monte Carlo 
realizations of IMBH populations in Milky Way-sized halos. 
Notice that unlike subhalo populations, IMBHs do not exhibit 
a significant anti-bias with respect to the DM. Rather, they 
are slightly biased toward being found near the halo center. 



rounding dark matter distribution in the vicinity of the 
IMBHs can survive tidal disruption to very small galac- 
tocentric distances. The final, cumulative radial distri- 
butions of unmerged IMBHs are shown in Fig. [21 They 
are very similar for scenarios A and B (though the nor- 
malization is different), and shows a behavior that scales 
as dN/dr ~ r~ 3 at large scales and tends toward a shal- 
lower slope on scales smaller than the scale radii of typical 
MW-size halos (see the following section). 



III. THE DENSITY ENHANCEMENT OF DARK 
MATTER AROUND IMBHS 

In each early-forming halo that hosts a seed black hole, 
when the black hole forms the DM distribution about the 
black hole inevitably reacts, adjusting to the new gravi- 
tational potential. This process has been studied exten- 
sively, particularly in the context of stellar cusps around 
massive black holes in clusters of stars or at the centers 
of galaxies (see e.g. |(52l Gondolo and Silk 

have applied this argument to the distribution of DM at 
the center of the Galaxy [7l| and introduced the term 
"spike" for the consequent enhancement in the DM den- 
sity around the central SMBH, in order to avoid confusion 
with DM "cusps" at the centers of halos in the cold dark 
matter model of structure formation. It was subsequently 
shown that dynamical processes like off-center formation 
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of the seed black hole, or major merger events, may lead 
to destruction or reduction of the spike [T^, uM ■ How- 
ever, steeply rising stellar cusps in the innermost regions 
of galaxies suggest that such processes were not effective, 
at least in the case of the Milky Way, or that the stellar 
cusps were re-generated via star formation or energy 
exchange between stars [75| . 

Recently, Bertone and Merritt studied the evolution 
of DM spikes including gravitational scattering off stars 
and the self-annihilation of DM particles [y, 0] , showing 
that the DM density in spikes is, indeed, substantially re- 
duced by these effects, but the enhancement of the anni- 
hilation signal is still significant with respect to ordinary 
DM cusps. 

In the present study, we are interested in "mini-spikes" 
surrounding IMBHs. Because we track the merger his- 
tory of each individual black hole, we can select precisely 
those black holes which never experienced mergers, to 
ensure that major mergers have not destroyed any cusp 
that existed around the original black hole. 

Furthermore, the models we explore predict from be- 
tween a few hundred to a few thousand black holes scat- 
tered throughout the Milky Way halo, and as the Milky 
Way has only 11 luminous companions within ~ 300 kpc 
|76|. we expect that the majority of the wandering black 
holes in our models reside in satellite halos with no sig- 
nificant stellar component. 

This implies that the effects of scattering off of 
stars should not significantly alter the DM distributions 
around the wandering IMBHs. The mini-spikes around 
unmerged, wandering, IMBHs are thus less sensitive to 
all of the dynamical processes that may have affected the 
spike at the Galactic center. 

We proceed now to evaluate the DM enhancements 
around IMBHs. As a first step, we need to specify the 
"initial" DM profile, that is, the DM distribution prior 
to black hole formation. Let the subscript "/" denote 
quantities at the time when the IMBH formed. The ini- 
tial DM profile of the mini-halo, before adiabatic growth, 
can be well approximated with a Navarro, Frenk, and 
White (NFW) profile |z3 

+ (2, 

The normalization constant po, and the scale radius r s , 
can be expressed in terms of the virial mass of the halo 
at the time when the IMBH formed M v [ r j, and the virial 
concentration parameter c V i r ,/ 

Cvir,/ Wir.f ) 

We recall that the virial mass is related to the virial ra- 
dius r v ir,/ by 

M vil -j = — [A vil (z f )p m (zf)] r%j , (4) 



while the function f(x) is, apart from constants, simply 
the volume integral of the NFW profile f(x) = ln(l + 
x)-\x/{l + x)}. 

In Eq. (0J, p m (zf) is the mean DM density at the red- 
shift of formation Zf, while A V j r (zj) is the virial over- 
density, for which we have adopted here the fitting form 
of Bryan and Norman [j^. At the redshifts of interest 
(z > 12) the universe is DM-dominated and the expan- 
sion rate and growth of perturbations are described by 
the standard relations for an Om = 1, "standard" CDM 
cosmology. In this case, A v i r (zf) ~ 18ir 2 ~ 178. For each 
black hole at redshift z — we extract from its merger 
tree the parameters M V i r ,/, c v i r ,f , Zf and use Eq. 
to calculate the initial DM profile before the formation 
of the black hole. Alternatively, we could have chosen 
the more recent parametrization proposed by Navarro et 
al. [73 (see also Refs. [8(], However, this profile 

implies modifications at scales smaller than those we are 
interested in, where the profile is anyway modified by the 
presence of the IMBH. 

We assume that the black holes form over a timescale 
long enough to guarantee adiabaticity, but short com- 
pared to the cosmological evolution of the host halo (in 
scenario B, both of these assumptions are built into the 
black hole formation model, see Section II B as well as 
our discussion below). 

Adiabaticity requires that the formation time of the 
black hole is much larger than the dynamical timescale at 
a distance from the black hole, where Th is the radius 
of the sphere of gravitational influence of the black hole, 
r/i ~ GA/bh/c 2 , and a is the velocity dispersion of DM 
particles at Th . In practice, we estimate by solving the 
implicit equation 

M (< r h ) = f " p(r)r 2 dr = 2 M hh . (5) 
Jo 

For a representative case in scenario B, with Mhh = 
1O 8 M and M vir ,/ = 1O 8 M , this gives r h /r s w 0.04. 
In scenario B, the black hole formation time is set by 
the timescale for viscous angular momentum loss and is 
limited by the evolutionary timescale of the first stars 
and the gravitational infall time across the gaseous disk, 
which is of order Myr (see Ref. 0| for a detailed dis- 
cussion of timescales). The relevant timescale for the 
mass build up of the IMBH is then t ev ~ 1 - 20 Myr. 
In scenario A, we follow Ref. 9] where the characteristic 
timescale for the growth of the black hole by accretion is 
taken to be of order 1 — 20 Myr for a plausible range of 
accretion efficiencies. 

The basics of adiabatic growth can be easily under- 
stood (e.g. Ref. 0), and in most cases the details can be 
worked out by taking into account the approximate con- 
servation of adiabatic invariants under a certain set of 
assumptions. If one starts from an initially uniform DM 
distribution, the final profile will be a mild mini-spike 
with density p sp oc (r/r/j) 3 / 2 (e.g. see [6^ and references 
therein). If one starts from a cuspy profile, such as the 
NFW profile of Eq. J2J), the new profile is essentially a 
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power-law, 



Psp(r) = p(r sp ) — 



(0) 



where the radius of the spike is r sp w 0.2rh l7p|._and 7 sp 
is related to the initial power-law index 7 by 



M, 8 

•y 7i] 



7s P 



9-27 
4-7 



(7) 



In the case of the profile of Eq. J2J) , this reduces to 7 sp = 
7/3. 

The DM annihilation flux in this case diverges at small 
radii. However, the very annihilations that we study 
here provide an upper limit to the DM number (and thus 
mass) density. In absence of other processes affecting the 
distribution of DM, the DM density obeys the equation 

n x (r,t) = -crv n 2 x (r,t) (8) 

where av is the annihilation cross-section times relative 
velocity (in the non-relativistic limit) and mx is the DM 
particle mass. The solution to the evolution equation is 



n x (r,t) 



l + n x (r,t f ) av (t-t f ) 



(9) 



which shows that efficent annihilations set an upper limit 
to the matter density of order m x /av(t — tf). We define 
Him as the radius where 

Psp(nim) = m x /av (t - t f ) = p\ im . (10) 

We therefore define an inner cut-off at a radius 



r cut = Max[4R Sch w,nim] 



(11) 



where i?s c hw is the Schwarzschild radius of the IMBH 
-Rschw = 2.95 kmMbh/ M©. For common values of the 
mass and cross section of the DM particle, ru m ~ 1CP 3 pc 
so that r cut = n im . 

Is the adiabatic growth of a central mass a good ap- 
proximation in our IMBH formation scenarios? We have 
already discussed the timescales involved, but the deriva- 
tion of the inner radius r\ ull provides us with the possi- 
bility of checking whether the size of the region where 
matter accretes, leading to the formation of a black hole, 
is actually smaller than the characteristic size of the DM 
spike, ri; m . In scenario A, this is not a problem, because 
in this case the black holes from from Pop. Ill stars and 
the spike is produced by the growth of the black hole, thus 
by processes occuring on scales of order i?schw << Him- 
In scenario B, the situation is different, because the mini- 
spike is produced by the flow of protogalactic material 
that lost its angular momentum by viscosity. 

Such an object may collapse directly to a black hole 
or it may form a short-lived, pressure-supported object 
[47ll63]| . However, in either case, the characteristic size of 
the massive object that forms is likely to be much smaller 



than rii m . We can make an order-of-magnitude estimate 
of the relative sizes as follows. Massive stars are believed 
to have a polytropic equation of state with n=3. In other 
words, the equation of state is described by 



P(r) = Kp{r) 1 



1 + 1/n, 



(12) 



and in this case n = 3 implies T — 4/3, as appropriate for 
a star supported by radiation pressure. It is possible to 
evaluate numerically the properties of a polytropic star 
in hydrostatic equilibrium, for n = 3 the approximate 
relation p c = 54. 2p between the central and the average 
density of the star holds (see e.g. |8^1- 

We infer that the typical scale for the radial extent of 
such an object should be 



■ 54 2 3M^ 
4irp c 



R* = ( 
io-V(n^)*(n=£ 



(13) 



This scale is clearly much smaller than the typical size of 
the spike r cut . 



IV. DARK MATTER ANNIHILATIONS IN 
MINI-SPIKES 

Dark matter particles are expected to have a non- 
negligible annihilation cross-section into Standard Model 
particles, in order to be kept in chemical equilibrium 
in the early Universe. It should be a weak or weaker- 
than-weak interaction in order to provide a relic den- 
sity which satisfies cosmological constraints (for recent 
reviews of DM candidates and detection techniques see 
e.g. Refs. 0, 0). Although it is difficult to make 
definitive statements on the nature of the DM parti- 
cles, it is commonly believed that a mass in the range 
m x ~ 100 — 1000 GeV would be a reasonable expec- 
tation in the most widely-discussed DM scenarios (e.g. 
minimal supersymmetry or scenarios with unified extra- 
dimensions). A naive estimate of the annihilation cross 
section, based on the observed relic abundance of DM, 
suggests that av ~ 10~ 26 cm 3 s _1 . This value can be 
more appropriately used as an upper limit to the an- 
nihilation cross-section, rather than an actual estimate, 
since processes like co-annihilations may significantly af- 
fect relic density yields (for more details see Refs. [l|, |2( 
and references therein). 

Instead of undertaking a detailed scan of the parame- 
ter space for different DM candidates, we limit ourselves 
here to estimates of the annihilation fluxes for two bench- 
mark models: an optimistic model, with m x = 100 GeV 
and av = 3 x 10~ 26 cm 3 s -1 , leading to large anni- 
hilation fluxes; and a model with m x = 1 TeV and 
av = 10~ 29 cm 3 s _1 , leading to more pessimistic predic- 
tions. We note that in both cases, the mini-spike profiles 
reach their maximum values at a radii ru m >> 4i?s c ii W , 
thus ri; m provides an estimate of the size of the region 
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where most of the annihilation radiation originates from. 
The case of annihilations from the DM spike at the center 
of the Galaxy has been extensively studied in the liter- 
ature in terms of neutrino ; gamma-ray, and synchrotron 
emission 0, H El S Hi IM . 

The flux of gamma-rays from a mini-spike around an 
IMBH can be expressed as 



$(E, D) 



1 ay I AN 
2mf Z) 2 " dE 



Psn( r ) r2dr 



dN 



Psp 



av r" 



dE 4 7s) 



6 D 2 



-27 sp +3 

UEi ) (14) 



where we assumed r sp >> r cut 



Inserting typical values 
of DM and spike parameters we get, for the case 7 = 1 
(7s P = 7/3), 



<5>(E,D) = $, 



dN 
3 dB 

D 
kpc 

' sp 

pe- 



ar 



3 cm 3 /s 



)( 



in v 



lOOGeV 

2 



10 2 GeVcm- 



10" 3 pc 



(15) 



with $o = 9 x 10 _10 cm _2 s _1 . It is useful here to empha- 
size the relative luminosities of IMBHs in the MW halo. 
In particular, consider the case of the relatively more lu- 
minous objects of scenario B. Using the fiducial values 
adopted in eq. I|15|). which are typical of scenario B, one 
can easily verify that the "luminosity" of a mini-spike 
(proportional to the volume integral of p 2 p ) is of the or- 
der of the gamma-ray luminosity of the entire Milky Way 
halo, a circumstance that has dramatic consequences for 
the prospects of indirect detection, as we describe in the 
following section. 

To estimate the flux, we need now to specify the 
gamma-ray spectrum per annihilation dN/dE, which de- 
pends on the nature of the DM particle. In most scenar- 
ios, direct annihilation in two photons is severely sup- 
pressed, but a continuum spectrum is expected from the 
decay of secondary neutral pions. 

In Fig. [3] we show the predicted gamma-ray spectra 
for different annihilation channels. For the bb channel 
we show two different curves, corresponding to differ- 
ent parametrization of the process of quark fragmenta- 
tion and subsequent decay of neutral pions, for 2 differ- 
ent mass scales of the DM particle. The first (FPS 04) 
corresponds to the parametrization in Ref. |87j, while 
the second (BSS 03) refers to the spectra presented in 
Ref. [g^ • The differences for different parametrizations 
and mass scales appear to be small, The third curve 
(BUB 97) corresponds to an analytic fit for the WW 
and ZZ channels, as discussed in Ref. jg^ , a channel 
leading to harder spectra with respect to the quark- 
antiquark channel. These channels often represent the 
most important annihilation channels for neutralino Dark 



E 



Z 
-d 



\ 



\ 



\: 



bb, rn^lTeV (FPS 04) 
bb. m,=100GeV (BSS 03) 
WW, ZZ (BBU 97) 
LKP (BBEG 05) 

0.1 

E/m y 



FIG. 3: Energy spectra of photons per annihilation for dif- 
ferent annihilation channels. The solid and dotted lines both 
correspond to the bb annihilation channel, the differences are 
due to different parametrizations of quark fragmentation and 
different DM particle mass scales. The solid line shows the 
parameterization of Ref. |87|| with m x = 1 TeV, while the 
dotted line shows that of Ref. with m x = 100 GeV. The 
short-dashed line corresponds to the spectra for annihilation 
through the WW and ZZ channels. In particular, we show 
the fit from Ref. [891 . Lastly, the long-dashed shows the spec- 
trum, summed over contributing channels, for annihilation of 
Kaluza-Klein DM from Ref. ||. 



Matter. In the case of Kaluza-Klein DM, other chan- 
nels become important. Following Ref. [j] we show in 
Fig- ED the total spectrum obtained by adding the contri- 
bution of different channels, weighted with the appropri- 
ate branching ratios. It is evident from the figure that 
in this case the spectrum is harder than the the quark 
or gauge bosons channels, due to contributions from in- 
ternal bremsstrahlung as well to decays of quarks and 
tau leptons. Internal bremsstrahlung is a general feature 
of scenarios where DM particles annihilate into pairs of 
charged fermions, which produces a sharp edge feature 
in the spectrum, dropping abruptly at a photon energy 
equal to the WIMP mass US Ell- 



in the next section, we will present our predictions us- 
ing the BSS 03 spectrum, with the caveat that different 
annihilation channels may lead to slightly different re- 
sults. The predictions for different annihilation channels 
can be easily obtained by plugging the appropriate spec- 
trum per annihilation into Eq. 1151 
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V. RESULTS 

In this section, we present our results on the detectabil- 
ity of annihilation radiation from the density enhance- 
ments about IMBHs. The number and properties of the 
IMBHs population are slightly different from one realiza- 
tion of Milky Way-sized halos to another, as described in 
Sec. II. To estimate the prospects of detection of IMBHs 
in the Milky Way, we thus need to average the results 
over all realizations. 

In Fig. we show the (average) integrated luminos- 
ity function of IMBHs in scenario B. We define the 
integrated luminosity function as the number of black 
holes producing a gamma-ray flux larger than $, as a 
function of <F The upper (lower) line corresponds to 
m x = 100 GeV, av = 3 x 10~ 26 crn^ 1 ( m x = 1 TeV, 
av = 10 -29 cm 3 s _1 ). In a practical sense, the plot shows 
the number of IMBHs that can be detected with experi- 
ments with point source sensitivity $ above 1 GeV. We 
show for comparison the point source sensitivity above 
1 GeV for EGRET and GLAST, corresponding roughly 
to the flux for a 5a detection of a high-latitude point- 
source in an observation time of 1 year [92] • The dashed 
region corresponds to the la scatter between different re- 
alizations of Milky Way-sized halos. This band includes 
the variation in spatial distributions of IMBHs from one 
halo to the next as well as the variation in the individual 
properties of each IMBH in each realization. 

Although one would naively expect that the fluxes 
scale with av/m x , we note that the DM profile itself de- 
pends on m x and av, more precisely on the ratio av/m x 
[see Eq. (|10|) ]. The maximum density is higher for the 
pessimistic case, and this partially compensates for the 
decrease in flux due to the prefactor av/m x . It is easy 
to see this in the case 7=1 from eq. H15|) as, by virtue 
of Ea. lTfjl r cut oc (<jv/m) -3 / 7 , and the final luminosity of 

the objects is thus proportional to ~ (av) 2 ^ 7 m x 9 ^ 7 . 

The number of detectable sources is very high, even in 
the pessimistic case, and either strong constraints on a 
combination of the astrophysics and particle physics of 
this scenario, or an actual detection, should be possible 
within the first year of operation of GLAST, which is 
expected to be launched in 2007. Depending on the spe- 
cific scenario, EGRET may have observed some of these 
IMBH mini-spikes, which would still account only for a 
small fraction of the unidentified sources. 

We show in Fig. [S] the integrated luminosity function 
of IMBHs in scenario A, for the same particle physics 
models shown in fig. |SJ The lines and error bars all have 
the same meaning as those in Figure 0] for scenario B. In 
this case of scenario A, mini-spikes are weaker, but the 
number of black holes is larger by roughly an order of 
magnitude, so that GLAST may still detect between a 
few tens and several hundred sources, whereas EGRET 
may have seen only a few or none. 

In figures ^ and [S] we have assumed that the main 
annihilation channel is bb. Although we have seen in the 
previous section that, depending on the nature of the 




lO-is 10" 1S 10-" 10-'° 10-» io-e 10 -7 10 -« 10-5 10- 4 10- 3 
*(E>lGeV) [cm" 2 s" 1 ] 



FIG. 4: IMBHs integrated luminosity function, i.e. num- 
ber of black holes producing a gamma-ray flux larger than a 
given flux, as a function of the flux, for our scenario B (i.e. 
for IMBHs with mass ~ 10 5 Mq). The upper (lower) line cor- 
responds to m x = 100 GeV, av = 3 x 10" 26 cm 3 s _1 (m x = 1 
TeV, av = 10~ 29 cm 3 s" 1 ). For each curve we also show the 
1-a scatter among different realizations of Milky Way-sized 
host DM halos. The figure can be interpreted as the num- 
ber of IMBHs that can be detected from experiments with 
point source sensitivity $ (above 1 GeV), as a function of <!>. 
We show for comparison the 5a point source sensitivity above 
1 GeV of EGRET and GLAST (1 year). 

DM particle, other channels may dominate and lead to 
different annihilation spectra. We see from these figures 
that the expected uncertainty, 0(1), would have a small 
influence on the number of objects that GLAST should 
be able to detect, certainly smaller than the uncertainties 
associated with m x and av and typically smaller than, 
or comparable to, the 1-a scatter between different Milky 
Way halo realizations. 

VI. DISCUSSION AND CONCLUSIONS 

We have studied the detectability of gamma-rays from 
DM annihilations in mini-spikes around IMBHs. The 
prospects of detection are summarized in figures 01 and 
where we show the number of IMBHs that can be de- 
tected from experiments with point source sensitivity $ 
(above 1 GeV), as a function of <&. We found that the 
prospects of detection with GLAST are so promising that 
a large number of sources may be detected within its 
first year of operation. With respect to the case of a 
spike at the Galactic center, searching for annihilation 
radiation from mini-spikes has the obvious disadvantage 
that IMBHs are smaller than the SMBH at the Galactic 
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10-'" lO" 13 10" 12 10-" 10- 10 10" 9 10" 8 lO" 7 10" 6 
*(E>lGeV) [cm" 2 s" 1 ] 



FIG. 5: IMBHs integrated luminosity function in scenario 
A (i.e. for IMBHs with mass ~ 10 2 M Q ). The upper (lower) 
line corresponds to m x — 100 GeV, av = 3 x 10~ 26 cm 3 s _1 
( m x = 1 TeV, av = 10 -29 cm 3 s _1 ). For each curve we also 
show the 1-cT scatter among realizations of Milky Way-sized 
halos. For the sake of comparison, we also show the point 
source sensitivity above 1 GeV for EGRET and GLAST. 



center, and the mini-spikes grow from less dense initial 
profiles. However, there are also several advantages. 

First of all, it is likely that the vast majority of mini- 
spikes around unmerged IMBHs in the outer Galactic 
halo are not affected by the dynamical processes that 
tend to destroy the central spike, or to decrease signifi- 
cantly the DM density near the black hole. For instance, 
they lack stellar cusps and are rarely affected by tidal 
interactions. Furthermore, the prospects of detectabil- 
ity appear very promising and certainly less problematic 
than, say, annihilations from the Galactic center. In fact, 
the gamma-ray background is strongest in the direction 
of the Galactic center and substantially reduced when 
observations are performed off the disk and in particular 
at high Galactic latitudes. There is even reason to be- 
lieve that the number density of IMBH may be enhanced 
at high Galactic latitudes |93|. 

Moreover, there are several known gamma-ray sources 
in the direction of the Galactic center, and the observa- 
tion of a unique source, even coincident with the Galac- 
tic center would not necessarily imply a DM annihilation 
origin. On the contrary, the detection of tens or more 
gamma-ray sources with identical spectra, in particular 
identical cut-offs at the DM particle mass, and not as- 
sociated with the Galactic disk or other luminous com- 
panions of the Milky Way, would provide a smoking-gun 
signature of DM annihilations. 

A natural place to search for IMBHs may be the known 



dwarf satellite galaxies about the Milky Way; however, 
the physics that govern the formation of these objects 
is still a topic of much debate and uncertainty, so such a 
search may be subject to the same drawbacks as searches 
for radiation near the center of the Galaxy. However, as 
we have already stressed, these IMBH formation scenar- 
ios have as a virtue that they predict that there may be 
hundreds of detectable objects within the Galactic halo, 
most of which would not be associated with the known 
population of dwarf satellite galaxies. An IMBH popula- 
tion similar to the one in the Milky Way halo should be 
present in Andromeda, given its similarity to the Galaxy 
in terms of size and mass. Moreover, the distance to the 
Andromeda IMBHs would be between « 400 - 1 000 kpc, 
amounting to a factor of only a few more in distance than 
to the black holes of the outer MW halo. Hence the de- 
tection of such IMBHs in Andromeda may be possible in 
optimistic scenarios and may serve to demonstrate the 
ubiquity of such phenomena. 

As a further implication, we found that the annihila- 
tion luminosity from any Milky Way-like halo may be 
dominated by annihilation around IMBHs in optimistic 
scenarios. As an example, we showed that in scenario 
B, the total luminosity of an individual mini-spike, in 
terms of annihilation radiation, may be comparable to 
the luminosity of the entire host halo. Therefore, such 
optimistic scenarios provide a significant "boost factor" 
for the gamma-ray background due to DM annihilations 
in halos at all redshifts (e.g., Ref. [94|, see also Ref. [9{| 
for a comparison between the prospects of indirect detec- 
tion from the Galactic center and the gamma-ray back- 
ground), as well as an enhancement of anti-matter fluxes. 
A detailed study of these alternative indirect searches re- 
quires a full analysis of the redshift evolution of IMBHs 
and spikes in halos of all masses in the former case, and 
of the propagation of anti-particles in the Galaxy in the 
latter, and it is thus beyond the scope of this paper. 

Interestingly, the prospects of indirect detection in this 
scenario do not depend strongly on the particle physics 
parameters. In fact, while e.g. in the case of annihi- 
lations from the Galactic center the annihilation flux is 
proportional to crv/m^., the flux from mini-spikes is lim- 
ited by the plateau in the number density due to DM 
annihilation itself For mini-spikes growing from 7 = I 
profiles, we have shown that the annihilation flux is in- 
stead proportional to (av) 2 ^ ra x ^ . 

Finally, we stress that the detection of these sources 
would shed new light on the origin of IMBHs and SMBHs. 
Mini-spikes appear to be ideal targets for large-field-of- 
view experiments such as GLAST. Another promising 
experiment with a large field-of-view could be the Alpha 
Magnetic Spectrometer (AMS-02) instrument, if prelim- 
inary estimates of its sensitivity to gamma-rays are con- 
firmed |9g|. Once the positions of the sources are de- 
termined, Air Cherenkov Telescopes such as CANGA- 
ROO II], HESS [13, MAGIC T3] and VERITAS 
may provide additional information, because of their bet- 
ter performance at higher energies and significantly bet- 
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ter angular resolution. The determination of a common 
cut-off in the spectra (possible only with ACTs for DM 
particles heavier than 300 GeV) will provide an estimate 
of the mass of the DM particle, while spectral features, 
such as annihilation lines or sharp edges, may provide 
important information on the nature of the DM particle. 
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